## read covariates (rasters)
fgdb_inputs <- file.path(pilot_path, "logReg_processedInputs.gdb")
## read template raster (environmental scores)
#env <- as.raster(arc.raster(arc.open(file.path(fgdb_inputs,"EnvExclusions_WECC_comb"))))
env <- raster(file.path(gdbtiffs,"EnvExclusions_WECC_comb.tif"))
names(env) <- "env"
####################################
## fix extents function
e <- extent(env)
fixExtent <- function(template, ext_template, in_ras){
## https://gis.stackexchange.com/questions/158159/snapping-raster-grids-in-r
in_ras = resample(in_ras, template, "bilinear")
in_ras <- crop(in_ras, ext_template)
if(sum(as.matrix(extent(in_ras))!=as.matrix(ext_template)) == 0){ # edited
in_ras_mask <- mask(in_ras, template) # You can't mask with extent, only with a Raster layer, RStack or RBrick
}else{
in_ras <- extend(in_ras, template)
in_ras_mask <- mask(in_ras, template)
}
return(in_ras_mask)
}
Suitable Wind Energy Project Areas: Western Region